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Abstract. The Brimstone Canary is a geographical highly variable specics, distributed throughout southern and eastern 
Africa. Here | present data on body measures and plumage characteristics of 476 skins from all parts of the species’ ran- 
ge (Fig. 1). In both univariate and multivariate analyses the wing and beak measures showed the most remarked diffe- 
rences between populations (Tab. 2, 5). Variation parallels climatic trends with individuals from hot and humid regions 
have shorter wings and smaller beaks (Tab. 7). In contrast to earlier studies plumage colour should not be used to distin- 
guish populations, because of high individual variation (Tab. 8). 


Only skins from South Africa (including Natal and Zululand) could be separated clearly (stronger beaks) from all other 
populations (shorter wings, less strong beaks). Clinal variation and high individual variation did not allow more diffe- 
rentiation in the northern part of the range. Therefore, the data are in agreement with the separation of the subspecies 
sulphurata (in South Africa) and sharpei (all other populations), while wilsoni and all other subspecies are not suppor- 


ted. 
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1. INFRODUCTION 


The traditional way to describe spatial variation of species 
is the recognition of subspecies. Subspecies are more or 
less clearly separated from other populations of the same 
species by differences in at least one trait and are distrib- 
uted in a more or less clearly defined subarea of the 
species’ range. This method gives, however, an incomplete 
picture of variation, because many traits vary in clines and 
the spatial characteristics of clines may differ from trait 
to trait (see GOULD & JOHNSTON 1972 and ZINK & REM- 
SEN 1986 for reviews). 

Hence, to distinguish a representative of a subspecies from 
an individual variant or a local form in a cline profound 
knowledge of individual and geographical variation 
within the species is necessary. Such knowledge was not 
available when most subspecies have been described and 
is stıll not available for many species. 


The Brimstone Canary Crithagra sulphurata ıs such a ge- 
ographical highly variable species. Within the Canaries 
five afrotropical species form the monophyletic genus 
Crithagra as suggested by morphological, ethological and 
molecular data (van den ELZEN 2000). C. sulphurata in- 
habits grasslands and savannas with scattered trees 
throughout southern and eastern Africa (Fig.1). Within this 


range and especially in its eastern part a high degree of 
morphological variation has led to the description of sev- 
en subspecies: 


— sulphurata (Linnaeus, 1766) (LINNAEUS 1766, Syst. Nat., 
ed. 12,1: 305), type from Cape of Good Hope, 

— sharpii (Neumann, 1900) (NEUMANN 1900, J. Ornithol. 
48: 287), t.t. Marangu, Kilimanjaro, 

- shelleyi (Neumann, 1903) (NEUMANN 1903, Ornithol. 
Mber. 11: 184), t.t. Kafuro, NW Tanzania, 

— frommi (Kothe, 1911) (KoTHE 1911, Ornithol. Mber. 19: 
71), t.t. Namanjera, Ufipa, SW Tanzania, 

— loveridgei (van Someren, 1921) (van SOMEREN 1921, 
Bull. Br. Ornithol. Club 41: 114), Lumbo, northern 
Mozambique, 

— wilsoni Roberts, 1936 (ROBERTS 1936, Ann. Transv. 
Mus. 18: 216), t.t. Kloof, Natal, and 

— languens Clancey, 1962 (CLANCEY 1962, Durban Mus. 
Nov. 6: 193), type from Manhica, Sul do Sava, south- 
ern Mozambique. 


Different authors have adopted controversial positions 
which of these types represent real subspecies, e.g.. FRY 
& KEITH (2004) in their recent review recognized only 
three of these (su/phurata, sharpei, wilsoni). RAND 
(1968) and CLANCEY (1972) reviewed available measure- 
ments and distribution data and concluded, in accordance, 
that most of the “subspecies” mentioned above are not 
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sharply separated from each other in the traits considered. 
Due to larger body size and stronger beaks South African 
nominate sulphurata can be distinguished from all other 
forms with relative ease. Therefore, RAND accepted only 
two subspecies: su/phurata and sharpei (all populations 
except South Africa). In an alternative view he divided 
sharpei on the basis of subtle differences in plumage 
colour and size into four subspecies: wilsoni (Natal. south- 
ern Mocambique), sharpei (Kenya), frommi (Tanzania, 
Zambia and Angola) and s/e//eyi (all other populations). 
In contrast, CLANCEY (1972) separated four subspecies in 
South Africa alone (su/phurata, wilsoni, languens and 
shelleyi). Both studies are, however, based on only a few 
traits and from some regions only a small number of in- 
dividuals have been included. 


In this study 1 present data on a larger variety of meas- 
ures and plumage characteristics of C. su/phurata skins 
from all parts of the species’ range. The aim of this study 
is to describe the extent and pattern of geographic varia- 
tion in this species. Which of the subspecies described ear- 
lier represent real taxonomic units and which are interme- 
diate stages arbitrarily picked out of a continuous cline of 
variation? In a second step 1 want to rclate morphologi- 
cal variation to possible causal factors such as ecological 
conditions (climate, vegetation, topography) and palaeo- 
graphic cvents. 


2. MATERIAL AND METHODS 


I measured a total of 485 muscum skins of Crithagra sul- 
phnrata, of which all 476 adults (253 males, 171 females 
and 52 unsexcd) were considered in the following analy- 
ses. Additionally, I studied 24 skins of the sister species 
(van den ELZEN 2000), Crithagra flaviventris, for compar- 
ison (15 males, 9 females). The skins are kept by Zoolo- 
eisches Forschungsmuseum Koenig (ZFMK), Bonn, Zo- 
ologisches Institut und Museum (ZMH), Hamburg, Zool- 
ogisk Museum (ZMUC), Copenhagen, Zoologische 
Staatssammlung (ZSM), Munich, Natural History Muse- 
um (BMNH), Bird Group, Tring, and Royal Museum of 
Central Africa (RMCA), Tervuren, respectively. 


Tarsus and beak measures (length, width and height in 
both cases) were taken from all specimens. Length of 
wing, tail and first to ninth primaries as well as Kipp’s dis- 
tance (wing tip to first secondary) and graduation of tail 
were measured in all except molting specimens. From the 
wing measures I calculated Wing Pointedness Index I and 
ll after MLIKOVSKy (1978) and HEDENSTROM (1989), re- 
spectively. Plumage colours (throat, breast, belly and back) 
were compared with standardized colour maps (Scandi- 
navian Colour Institute 2001). For statistical analysis I de- 
scribed standard colours by four variables: proportion of 


black, of colour, green (negatively correlated to propor- 
tion of yellow due to construction of the colour system), 
and red (C. flaviventris only). Not all traits could be tak- 
en for all individuals because of skin condition. Therefore, 
sample sizes differ from analysis to analysis. 


Skins with geographical origin close to each other have 
been pooled for analysis, resulting in 15 ‘populations’ (Fig. 
1). Individuals which locality information were lacking, 
unclear or doubtful (e.g. Collection Meinertzhagen; cf 
KNOX 1993) were included in the analyses, but considered 
as belonging to none of the populations. To get larger sam- 
ple sizes some of the populations were pooled to larger 
‘population groups’ (A-G) in some analyses (Fig. 1). 
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Fig. 1. Map of collection localities of skins (crosses - single 
skins; numbers - series with number of collected individuals). 
Population groups” are marked 


+9 Gk 


Ellipses define “populations”. 
by letters A to G. The present breeding distribution of C. sulp- 
hurata (simplified after FRY & KEITH 2004) is shaded. 


Predominant vegetation structure was revealed for each 
skin locality using the vegetation map of EC-JRC (2003, 
map 2) summarizing all categories under “dominant tree 
layer’ as ‘forest’ and all categories under ‘dominant agri- 
culture’ and ‘dominant shrub or grass layer’ as ‘open habi- 
tat’. Climatic data of the closest weather station were tak- 
en from LEBEDEV (1970). Because of the scarcity of weath- 
er stations, up to 63 skins gathered on one weather sta- 
tion and, in these cases, I used the mean for each trait. The 
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following data were considered: yearly mean of mean, 
maximum, minimum and range of air temperature per day, 
yearly sum of precipitation, daily maximum of precipita- 
tion. Altitudinal effects could not be taken into account 
because for most skins altitudinal information were not 
available with sufficient accuracy. 


Non-parametric Factor Analysis were performed using the 
MS Excel Add-In XLSTAT 7.0 (ADDINSOFT 2003). For all 
other calculations I used the software package SPSS for 
Windows 11.0 (SPSS 2001). 


3. RESULTS 
3.1. Data preparation 


Some variables were transformed to fit a normal distri- 
bution: beak height (reciprocal square root transforma- 
tion), first primary length (reciprocal), graduation of tail 
(natural logarithmic) and the tarsus measures (cube). All 
other metric variables did not differ significantly from nor- 
mal (Kolmogorov-Smirnov One-sample-test, p>0.05). 


The plumage colour variables have been subjected to a 
non-parametric Factor Analysis based on Spearman’s rank 
correlation. After Varimax rotation the resulting three fac- 
tors with eigenvalues larger than | explain 88 % of the 
variance (Tab. 1). The first factor (hereafter called ‘front 
colour’ or ‘C1°) is highly associated with belly and throat 


Table 1. Correlations between plumage colour variables and 
factor scores from a Factor Analysis based on a Spearman’s rank 
correlation matrix. Correlations > 0.6 printed bold. 





Factor Cl 2 C3 

belly: % black -0.367 0.226 _ 0,027 
belly: % colour 0.751 -0.053 -0.134 
belly: % green -0.740 0.045 0.175 
breast: % black -0.056 0.248 0.943 
breast: % colour 0.336 -0.223 -0.747 
breast: % green -0.225 0.010 0.620 
back: % black -0.147 0.717 0332 
back: % colour 0.243 -0.951 -0.188 
back: % green -0.141 -0.307 0.026 
throat: % black -0.567 OTO "OTIS 
throat: % colour 0.723 0.002 -0.145 
throat: % green -0.776 0.015 0.248 
eigenvalue 4.4 1.6 1.1 

% variance explained Son 223 27.8 


colour; higher valucs represent individuals with morc yel- 
lowish and less greenish front plumage. C2 (“back colour’) 
describes back colour, higher values correspond with morc 
black and less colour components in the back plumage. 
The colour of the breast is highly correlated with C3 
(“breast colour’). Individuals with high values of C3 have 
darker and more greenish breasts while low C3 values in- 
dicate yellowish breasts. For further analysis I used the 
factor scores of C1-C3 for each skin. 


3.2. Morphological differences between C. sulpliurata 
and C. flaviventris 


At first I looked at differences in morphological charac- 
teristics between C. su/phurata and its sister species, C. 


flaviventris. A Discriminant Function Analysis (DFA) re- 


duces the 21 morphological variables to a single canoni- 
cal variable with an eigenvalue of 1.13 and a canonical 
correlation of 73 %. Two beak measures and two wing 
measures are included in the analysis (standardized canon- 
ical coefficients in parenthesis): beak length (1.533), beak 
height (0.696), Ist primary length (0.463) and distance 
from wing tip to secondaries (0.343). Valucs for C. fla- 
viventris are lower, thus, skins of this species are small- 
er in beak and wing dimensions. 95 % of all skins have 
been classified correctly by the DFA. An additional trait 
for differentiation is a red component in the front 
plumage of C. flaviventris: 13 of 21skins had a slight tinge 
of red, but none of the C. sulphurata skins. 


3.3. Differentiation of C. sulphurata populations 


Most traits differ significantly between populations (Tab. 
2). Wing and beak measures show high levels of differ- 
entiation (high significance values) in both males and fe- 
males. For a detailed analysis I tested for each trait if there 
is a difference between population groups. Again, the most 
pronounced differences are in the wing and beak meas- 
ures (Tab. 3). Of 21 population group pairings 16 (males) 
and 11 (females) differ significantly in at least one wing 
measure, respectively, and 12 (males) and 6 (females) in 
at least one beak measure, respectively. The populations 
of Kenya (group G) and South Africa (A) are separated 
by differences in most traits from all other population 
groups, but not from each other. In the central part of the 
species’ range populations vary hardly. There is a tenden- 
cy for females to differ in fewer traits than males across 
the same population group pairings. 


Many of the studied traits correlate with each other (e.g. 
wing length and primary length). Intercorrelations may oc- 
cur between all variables. Therefore, I performed a Prin- 
cipal Component Analysis (PCA) on all 24 original meas- 
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Table 2. Univariate differences in C. sulphurata traits between populations (one-way ANOVA). *significant differences on the 
5%-level (corrected for 24 tests: P<0.05/24=0.0021). 


males females 

df F Sig. df E Sig 
WINQ MEASUFCS 
wing length 13215 19.002 0.000* 13/146 10.548 0.000* 
Ist primary length 13/192 8.443 0.000* 11/127 4.040 0.000* 
2nd primary length 13/192 7.286 0.000* 114128 4.660 0.000* 
3rd primary length 132197 7.083 0.000* 112129 6.708 0.000* 
4th primary length 13/189 Tease 0.000* 11/124 4.818 0.000* 
Sth primary length 13/195 11:383 0.000* 11/129 6.202 0.000* 
6th primary length 13/194 15.328 0.000* a) E72] 0.000* 
7th primary length 13/195 19.416 0.000* 11.182 11.504 0.000* 
Sth primary length 13/188 17.909 0.000* 11/131 10.368 0.000* 
9th primary length 13/191 1272] 0.000* 11/130 7.282 0.000* 
Kipp 13/213 2.797 0.001 * 13/142 1.433 lei 
wing index I 9/61 0.534 0.844 9/39 1.242 0.299 
wing index I] 9/61 0.494 0.873 9/39 1.082 0397 
fail measures 
length 13/216 7.475 0.000* 13/145 5.812 0.000* 
graduation 13/214 1.230 0.299 13/143 12025 0.431 
beak measures 
length 137217 9.009 0.000* 13/147 2020 0.000* 
height 12/104 11.610 0.000* 13/73 2 0.000* 
width 13/218 20.355 0.000* 12/146 IZOLO 0.000* 
tarsus Measures 
length 13/216 4.631 0.000* 13/147 3.1906 0.000* 
width 13/215 2,9706 0.002 13/147 2.681 0.002 
height 13215 2.650 0.002 13/147 3,032 0.000* 
Plumage colour 
front colour 13/218 3.109 0.000* 13/147 1.487 0.129 
back colour 13/218 4.592 0.000* 13/147 2.047 0.021 
breast colour 13/218 4.067 0.000* 13/147 3.002 0.000* 


ures (21 morphometric and 3 colour traits). It reduces the dividuals with stronger beaks (beak height transformed as 
data set to five uncorrelated factors that together explain lst primary). PC3 and PC4 represent wing shape and 
65 % of the variance (Tab. 4). The correlations between breast colour, respectively. The fifth factor correlates neg- 
the first factor (PC1) and length of wing, tail and all pri- atively with front plumage and positively with back 
marics are all relatively large and positive (lst primary ap- plumage. Individuals with high PCS scores are, therefore, 
pears inverse due to transformation). It can be taken as a darker on their back and more greenish on their front, 
multivariate measure of size; individuals with high PCI while individuals with low PCS scores are more greenish 
values are larger in size. The second factor summarizes on their back and more yellowish on their front. 

the three beak measures: higher PC2 values indicate ın- 
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Table 3. Differentiation of traits (raw measures) between the population groups A to G. Shown arc traits with significant diffe- 
rences on the 5%-level (Mann-Whitney-U-test, significance level corrected for 24 tests as in Tab. 2). The diagonal separates ma- 
le (right) and female (left) analyses. 

Abbreviations: | — length, w — width, h — height; 1-9 first to ninth primary length. B, C, G: proportion of black, colour and green 
of breast, throat, belly and back plumage, respectively. 
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Table 4. Correlations between the 24 original measures and principal component scores from a PCA (using Varimax rotation with 
Kaiser normalization) of the correlation matrix of 476 C. sulphurata skins. Correlations >0.600 are printed bold. 


PC? 


Variable PCI PC3 PC4 PCS 
wing length 0.864 0331 0137 0.105 0.013 
Ist primary length -0.726 -0.171 0.199 -0.108 -0.156 
2nd primary length 0.780 0.124 -0.274 0.087 0.175 
3rd prımary length 0.800 0.120 -0.292 0.029 0.120 
4th primary length 0.822 0.077 -0.298 -0.013 0.059 
Sth primary length 0.858 0.119 -0.236 0.039 0.011 
6th primary length 0.891 0.208 -0.012 0.103 0.022 
7th prımary length 0.919 0215 0.091 0.100 0.040 
8th primary length 0.905 0.224 0.141 0.092 0.050 
9th primary length 0.844 0.203 0.137 0.100 0.102 
Kipp distance 0.366 0.198 0.234 0.036 -0.208 
wing index I -0.111 -0.029 0.946 0.037 0.041 
wing index H -0.110 -0.012 0.943 0.038 0.007 
tail length 0.636 0.273 0.081 0.125 0.105 
tail graduation 0.012 0.182 -0.035 0.461 -0.082 
tarsus length 0.286 0.048 0.052 0.625 0.189 
tarsus width 0.094 0.369 -0.015 0.106 0.296 
tarsus height 0.181 0.294 0.004 0.209 oo 
beak length 0.400 0.697 -0.010 0.148 -0.009 
bcak height -0.218 -0.823 0.033 -0.020 -0.023 
bcak width 0.382 0.729 -0.004 0.078 0.096 
tront colour 0.032 -0.153 0.027 0.133 -0.661 
back colour 0.170 -0.055 0.050 0.009 0.683 
breast colour 0.036 -0.011 0.071 0.810 -0.081 
cıgenvalucs 8.172 2.463 Ze 1.452 1.296 
% variance explained 34.0 10.3 9.5 6.1 5.4 


Table 5. Standardized Canonical Discriminant Function Coefficients from a DFA of the five PCA scores of 443 C. sulphurata 


skins from the population groups A to G. 


PCA-factor DFI DF2 DF3 DF4 
PCI 6.824 -0.633 -0.272 -0.007 
PC? 0.816 0.612 -0.32 -0.184 
PC4 0.52 -0.026 0.779 -0.455 
PES 0.516 0.178 0.365 0.809 
eigenvalues L517 0.244 0.136 0.015 
% variance explained 793 [27 7.1 0.8 
canonical correlation 0.443 0.346 


0.776 


0.122 
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Table 6. Proportion (%) of skins classified correctly as belon- 
ging to population group A-G by the DPA. N - sample size. 


To determine which of these factors have the highest sc- 
lectivity for the population groups, | performed a stepwise 
Discriminant Function Analysis (DFA) on the PCA scores 
(Tab. 5). PC3 (wing shapc) does not contribute to differ- 


B D E F G entiation between groups and is removed by the analysis. 
N 29 79 0 0 0 3 a i The other ut PCA scores are transformed to four DFA 
factors. 65 % of skins are classified correctly by this analy- 
B 21 0 0 0 0 0 100 0 sis. There are, however, marked differences between the 
C 36 0 0 0 0 8 89 3 population groups in the proportion of correct classifica- 
D 14 0 0 0 0 >] 7 0 tion (Tab. 6). Population group A (South Africa) is clear- 
ly divided off from all other groups. Population group G 
E 72 3 0 0 0 40) 53 4 T e f 0 RR 
| (Kenya) is separated less, but more than 50 % of skins are 
eel l 0 l 0 3 9 1 classified correctly. Most overlapping of G-skins is with 
Guo 34 9 0 0 0 18 2} 53 group A. Nearly all other skins from groups B, C, D, E 
- -- and F are classified as belonging to population group E 
(Congo, SW Tanzania) or F (Kivu/Uganda/Rwanda). Sep- 
arating males and females in a single DFA each gives sim- 
ilar results and is, therefore, not presented here. 
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Fig. 2. Crithagra sulphurata skins in the space defined by principal components | and 2 (see Table 4). Solid symbols - all skins; 


filled symbols - skins from population group A to G, respectively. 
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The groups B to F of eastern Africa which appear homog- 
enous in the above population group analysis may be dif- 
ferentiated on a lower level. Therefore, I calculated a DFA 
at population level on all skins except from population 
group A and G. Again, most skins are classified as belong- 
ing to populations of group E or F. 


3.4. Pattern of geographic variation 


The large range of classification accuracy by DFA (Tab. 
6) indicates that some population groups are more sepa- 
rated from their neighbours than others. Scores of all 
groups overlap to some extend with the characteristics of 
other populations (Fig. 2). The pattern of variation 1s much 
more complex than described by clear-cut population 
groups. As an example, in Fig. 3 / Fig. 4 wing length and 
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Variation of wing length of C. sulphurata. Shown are 


line schematizes the coastline of southern and eastern Africa. 


beak width are shown against the geographical origin of 
the skins. 


The population group with lowest overlap to other groups 
is South Africa (including Natal; group A). South African 
C. sulphurata are characterized by strong beaks and large 
to intermediate body sizes (Fig. 2). Within population 
group A both wing and beak measures are relatively ho- 
mogeneous (Fig. 3, Fig. 4). Skins from Natal have a lit- 
tle bit lower scores ın both principal components, thus dis- 
tinguished by little weaker beaks and little shorter 
wings/tails. The differences to the next population in the 
north (group C) are, however, much more remarkable than 
the differences between the South African and Natal pop- 
ulation (Figs 2—4). 
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In relation to all other groups Kenyan C. sulphnrata (group 
G) look similar to skins from the other end of the species’ 
range, group A. Beaks of skins from Kenya are interme- 
diate between South African birds and all other groups, 
while body size is larger in Kenya than in all other groups 
(Fig. 2). Overlap to neighbouring populations is, howev- 
er, considerable. As shown by Fig. 3 and Fig. 4 the pop- 
ulation of Kenya is less homogenous in both beak and 
wing measures. Skins from northern Kenya have as long 
wings as skins from southern Kenya, but in respect to beak 
size they resemble birds from the eastern and central part 
of the species’ range. 


Outside Kenya and South Africa no remarkable morpho- 
metric differentiation between population groups exist. 
There is a tendency, however, for skins from Mocambique 
(B) and Angola (D) to be smallest in both principal com- 
ponents (Fig. 2). Wing length increases from the coasts 
to central Africa and could be described as a quadratic 
function of geographic longitude (WING=-1.23E- 
10*LONG? + 6.87E-05*LONG + 66.25; R?=0.06; ANO- 
VA df=2/374 F=11.1 P=0.000). Beak width variation is cli- 
nal, too, and fits a linear function with highest values at 
the eastern coast (BEAK=-1.51E-06*LONG+6.82; 
R2=0.03; ANOVA df=1/376 F=11.2 P=0.001). Apart from 
this clear trend, for both measures areas with low scores 
alternate with high score areas, however (Fig. 3, Fig. 4). 


The clinal longitudinal trend is visible in the PC1/PC2 plot 
as a shift along the axis B (Mozambique) to C 
(Malawi/Zimbabwe) to E (Congo/Tanzania; Fig. 2). Due 
to high intra-group variability both population group F 
(Kivu/Rwanda/Uganda) and D (Angola) overlap with 
nearly all other central population groups. Skins from 
these two groups are, therefore, indistinguishable from the 
other population groups except A and G. Hence, the DFA 
could not dissolve the population groups B to F (Tab. 6) 
and classified all skins from the central part of the species’ 
range as belonging to group F. 


3.5. Impact of climate and vegetation 


To reveal associations between variation of morphologi- 
cal traits of C. sulphurata and climatic conditions | per- 
formed a PCA with the six climatic variables mentioned 
above and the five principal components of morphologi- 
cal measures (cf. Tab. 4). The resulting two factors explain 
45 % of total variance (Tab. 7). The first principal com- 
ponent (cPC1) shows that body size (PC1) decreases along 
a gradient of increasing temperature. Large C. sulphura- 
ta can be found under relatively moderate temperature 
conditions. Beak dimensions decrease as sum of precip- 
itation increase. Therefore, strong beaks are associated 
with relatively dry climate. 


Table 7. Correlations betwcen the 6 climatic variables, the 5 
principal components (cf. Table 4; mean for each weather sta- 
tion) and principal component scores from a PCA (using Vari- 
max rotation with Kaiser normalization) of the correlation ma- 
trix of 68 weather stations. Correlations >0.6 are printed bold. 


variable cPC 1 cPC2 
mean of temperature 0.898 -0.240 
maximum of temperature 0.850 -0.146 
minimum of tempcrature 0.916 -0.089 
range of temperature -0.262 -0.057 
sum of precipitation -0.041 -0.801 
maximum of precipitation -0.468 -0.064 
PC] -0.641 0.443 
PEZ 0.149 0.785 
PES -0.156 0.016 
PC4 -0.089 0.327 
PC5 -0.238 0.329 
eigenvalues 3.4 LS 
% variance explained 31.0 13.9 


Table 8. Mean variance of principal components (scc Tab. 4) 
at eight localities (“local”; between Sand 15 skins, total N=59) 
and in a random sample of the respective populations (“popu- 
lation”; same sample sizes). “Significant”: number of the eight 
local-population pairings with significant differences in varian- 
ce (Levene s test, p=0.05). 


Component | Local Population Significant 
PC1 (body size) 0.59 0.93 2 
PC? (beak size) 0.48 0.79 | 
PC3 (wing shape) 1.03 0.53 l 
PC4 (breast plumage) 0.77 1.19 l 
PC5 (front/back plumage) 0.80 1.44 l 


Skins from forested areas have lower scores of the sizc- 
related first principal component (thus, are smaller in size: 
mean -0.12) than skins from open landscape (mean 0.07; 
t-test t=2.06 df=416 P=0.040; P all other PC >0.05). This 
difference is more pronounced in Kenyan skins (popula- 
tion group G), where it is highly significant (t=3.76 df/=32 
P=0.001). 
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3.6. Individual variation 


From seven localities J have series of eight or more skins 
that had been sampled in the same season (within three 
month). I compared variability of these samples with vari- 
ability in a random sample of same sample size taken from 
the same population. In principle variances of local sam- 
ples are slightly smaller than population variances (Tab. 
8), but this difference 1s significant in only one or two of 
eight pairings. Mean variances of the eight samples are 
smaller ın principal components | and 2 (body and beak 
sizc) than in the other variables (plumage colour). This is 
true for both local and population level samples. 


4. DISCUSSION 


The morphological variability of C. sulphurata is consid- 
erable. Differences in body dimensions and plumage 
colour occur not only between populations but also indi- 
vidually at a single locality (Tab. 8). Therefore, ıf geo- 
graphical variation is studicd attention has to be paid to 
individual variation, too. 


In both univariate and multivariate analyses wing and beak 
measures show the most striking differences between pop- 
ulations (Tab. 2, 5). Wing and beak measures also cxhib- 
it lowest individual variation at one locality (Tab. 8). 
Therefore, these traits are suitable to separate populations 
and subspecies. The revisions of RAND (1968) and 
CLANCEY (1972) werc based on wing and beak measures, 
too. Additionally, both authors considered slight differ- 
ences in plumage colouration (a little more greenish/ycl- 
lowish etc.) to recognize subspecics. My study of a larg- 
er sample of skins shows, however, that plumage colour 
varies enormously at one site (Tab. 8) and can not be used 
to distinguish populations. 


Morphological variation of C. sulphurata parallels climat- 
ic trends: individuals from hot and humid regions have 
shorter wings and smaller beaks (Tab. 7). The pattern of 
variation in wing and beak size (Fig. 3, Fig. 4) could be 
explained by altitudinal effects to a large extend. Similar 
climatic conditions may be responsible for similar trait di- 
mensions of the populations of Kenya and South Africa, 
resident at the opposite borders of the species’ range. Dif- 
ferences in body size ın relation to predominate vegeta- 
tion cover are probably a side-effect of climate, too, be- 
cause also in rain forest areas C. sulphurata settles in the 
scattered open habitats and not in forests (FRY & KEITH 
2004). Especially in Kenya, where C. su/phurata inhab- 
its a broad range of altitudes, the relation between body 
size on one hand and climate and vegetation on the oth- 
er hand is strong. In this population separation of sites 
close to each other but differing in altitude is obvious. If 


wing length is taken as a measure of body size, this cli- 
matic trend parallels Bergmann’s rule, which suggest an 
increase of body size from cold to warm climates (ZINK 
& REMSEN 1986). 


There 1s little evidence for geographical isolation of pop- 
ulations due to climatic changes in the past 20.000 years. 
Palaeovegetation records certify distinct changes ın veg- 
etation cover during this period (ADAMS 2004). Habitats 
occupied by C. sulphurata today (grasslands, scrub, sa- 
vannas) have been driven back due to moister conditions 
during the holocene (about 10.000 years ago), resulting 
probably in a much smaller range of C. sulphurata at that 
time. However, open landscapes ranged continuously from 
southern A frica to Kenya, at the most interrupted by a nar- 
row forest belt at 20-25° S. 


Geographic variation reflects adaptation to different en- 
vironmental conditions within a species’ range. Therefore, 
it could be seen as a model of evolutionary steps of a sin- 
gle population over time in a changing environment 
(GOULD & JOHNSTON 1972). Differences in wing shape and 
size between populations often parallels migratory behav- 
tour (LEISLER & WINKLER 1985). FRY & KEITH (2004) de- 
scribe C. sulphurata as resident throughout its range. Short 
distance (altitudinal?) movements may occur in some pop- 
ulations, but at the most occasionally (MACKWORTH-PREAD 
& GRANT 1963, HARRISON et al. 1997). Therefore, mor- 
phological differences more probably reflect differences 
in habitat structure and usc since body size and shape af- 
fect manocuvrability (ZINK & REMSEN 1986). As LEISLER 
& WINKLER (1985) have shown ın their comparative study, 
morphological characters are strongly correlated with for- 
aging techniques. The variability of beak dimensions in 
this study suggests differences in foraging strategies or 
food composition of habitats within the range of C. sul- 
phurata. In all parts of the range C. su/phurata settles in 
open habitats with scattered trees, food mainly consists of 
seeds and small fruits (MACKWORTH-PREAD & GRANT 
1963, HARRISON ct al. 1997, FRY & KEITH 2004). How- 
ever, detailed information on geographic variation of habi- 
tat and food preferences is lacking. More field studies of 
the species are necessary to reveal this relationship. 


Single traits can be linked due to developmental con- 
straints. The PCA shows that the different measures of 
wing and tail size do not differ independently, but repre- 
sent a single principal component (Tab. 4). This corre- 
sponds with many other studies, where the first, highly 
correlated factor could be interpreted as a size factor 
(GOULD & JOHNSTON 1972). The fact, that beak dimen- 
sions, Wing shape and plumage colour each form separate 
factors, points to habitat or food related adaptations inde- 
pendent of size. 
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Only the skins from South Africa can be separated clear- 
ly from all other populations with the data presented here. 
The clinal nature of variation, the large overlap of trait 
measures and the large individual variation do not allow 
the recognition of subspecies in the remaining, northern 
part of the species’ range. 


C. sulphurata of the South African nominate subspecies 
sulphurata are characterized mainly by their strong beaks 
(Fig. 2). The range of this subspecies includes South Africa 
from the Cape Province to Natal and Zululand. In contrast, 
RAND (1968) recognized the population of the Cape 
Province as a separate subspecies, too, but classified the 
populations of north-eastern South Africa (Natal/Zululand) 
as belonging to another subspecies which includes all oth- 
er populations in the north. Already MEEs (1970) point- 
ed to the fact, that the populations of the Cape Province 
and Natal are different from each other in plumage colour 
at the most, but not in other measures (see also RAND’s 
own measurements). He therefore argued for combining 
both populations into one subspecies. The statistical analy- 
sis of skins presented here allows no separation between 
individuals from the Cape Province and Natal/Zululand, 
but separates between these and the populations in the 
north. These are isolated geographically, too, because C. 
sulphurata is lacking along the Limpopo River (HARRI- 
SON et al. 1997). 


Individuals distributed in the northern part (sharpei) are 
— despite marked individual variation — smaller (shorter 
wings) and have less strong beaks than their South African 
congeners. Both traits vary outside South Africa in a com- 
plicated pattern without abrupt changes. Variation in both 
wing length and beak size have a clinal and an altitudi- 
nal compound. The separation of subspecies 1s arbitrari- 
ly and would not contribute to a better understanding of 
geographical variation in this case. At best the population 
of the Kenyan highlands could be recognized by relative- 
ly strong beaks, but the transition to neighbouring popu- 
lations is fluid. A possible differentiation in traits not stud- 
ied here (e.g., genetical or ethological) have to be reserved 
for future studies. 


In conclusion, the data presented here are in agreement 
with the viewpoint of FRY & KEITH (2004) in respect to 
the subspecies sulphurata and sharpei, while the recog- 
nition of a subspecies wi/soni is not supported by this da- 
ta set. A further differentiation of populations in the north- 
ern part of the range as done by other authors (see above) 
could not be confirmed in this study. 
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